About matter and dark-energy domination eras in i?" gravity or lack thereof 
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We provide further numerical evidence which shows that R" models in f{R) metric gravity 
whether produces a late time acceleration in the Universe or a matter domination era (usually 
a transient one) but not both. Our results confirm the findings of Amendola et aZ. [ij-Qli but using 
a different approach that avoids the mapping to scalar-tensor theories of gravity, and therefore, 
dispense us from any discussion or debate about frames (Einstein vs Jordan) which are endemic in 
this subject. This class of models has been used extensively in the literature as an alternative to the 
dark energy, but should be considered ruled out for being inconsistent with observations. Finally, 
we discuss a caveat in the analysis by Faraoni [Jj , which was used to further constrain these models 
by using a chameleon mechanism. 

PACS numbers: 04.50.Kd, 95.36.+X 



I. INTRODUCTION 

f{R) theories of gravity are perhaps the most straight- 
forward modification of general relativity (GR), provid- 
ing an extra geometric component which in some partic- 
ular cases is capable of generating the accelerated expan- 
sion of the Universe manifested in supernovae la [5|. A 
large amount of literature has been accumulated in the 
past ten years about this kind of alternative theories of 
gravity and is beyond the scope of the present article to 
make justice to this vast subject (see Refs. ^ for a thor- 
ough review). Although some specific f{R) models have 
shown to be consistent with certain astronomical obser- 
vations, within the Solar System and also cosmological, 
not every model has the same success, for instance the 
/(i?) = Ai?„(i?/i?„)" model, simply referred in the litera- 
ture as to i?". Recently, Amendola et al. [ll,[l| performed 
a detailed analysis on the cosmological viability of several 
classes of f{R) models, including i?". Using a dynami- 
cal system approach, they concluded that for this latter 
the usual matter era that precedes the accelerated phase 
with an scale factor a{t) ^ t^/^ is generically replaced by 
an non standard era with a{t) - (c.f. Ref. for 
a complementary analysis), and in the cases where it is 
possible to achieve a usual matter domination epoch the 
accelerated expansion is not possible. In any instance, 
the conclusion was that such a model is simply unable to 
reproduce the observed features of our Universe without 
the addition of some form of dark energy. 

These results have been, however, the object of a de- 
bate concerning two issues: 1) the frames (Einstein vs 
Jordan) used in the scalar-tensor (ST) approach to ana- 



lyze the i?" and other models and 2) the analysis 

of the phase space [13, [HI ■ 
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Since the i?" model has been and keeps being consid- 
ered in the literature (see a complete list of references 
in [3|) it is important to settle this question with an in- 
dependent method an beyond any reasonable doubt. 

In this brief report we reanalyze the cosmological case 
of the i?" model using a different method and spanning 
a wide range of n. Our technique does not involve what 
is usually called the scalar-tensor approach (ST) where 
a scalar field (j) = fa is defined in order to map f{R) 
theories to a Brans-Dicke like theory with uj — and a 
potential. Instead, we promote the Ricci scalar itself as 
an independent degree of freedom [l2l . IT3| and in this way 
we circumvent the potential drawbacks associated with 
the ST approach (e.g. multivalued scalar-field poten- 
tials) , and in addition avoid the long standing issue about 
frames (Jordan vs Einstein) which plagues not only the 
ST method, but also the analysis of scalar-tensor theories 
themselves, and which gave rise precisely to the unneces- 
sary debate mentioned above about the cosmological vi- 
ability of i?" and other class of f{R) theories. As we will 
show, our approach leads to a rather "friendly" system 
of equations which are much more simple to treat than 
other systems found in the literature and that can be eas- 
ily solved numerically. We had used this method before 
in the analysis of compact objects fl^l and more recently 
in cosmology using different f{R) models [H, [l3l- For 
the cosmological analysis at hand, we shall consider the 
same tools developed in fl3| and adapt them to the case 
R'\ 

Our analysis supports the general conclusions of [l|, [^ 
and (3 (although we do no commit ourselves in assessing 
the soundness of their phase-space analysis) providing 
a second, independent, strong and unambiguous piece 
of evidence showing that the specific i?" model is not 
cosmologically viable. In the next section we discuss in 



detail our findings that lead to such conclusion, and we 
also argue that the analysis put forward by Faraoni Q 
to constrain these kind of model in the light of the Solar 
System tests using a chameleon mechanism, is ill founded 
and requires a deeper review. 



II. /(i?) THEORIES 



The action in f{R) gravity is given by: 

rfiR) 



S[9ab, ■0] 



2k 



gd X + Smuttigabilp] , (1) 



where k = SttGq (we use units where c = 1), f{R) is a 
sufficiently differentiable but otherwise a priori arbitrary 
function of the Ricci scalar R and represents schemat- 
ically the matter fields. The field equation obtained from 
Eq. (m is: 



fliRab 



f9ab - (VaVb - gab^) Ir ^ K.Tab , (2) 



where fa indicates dnf, □ — g'^'^^a^b is the covariant 
D'Alambertian and Tab is the energy-momentum tensor 
of matter associated with the xjj fields. From Eq. ^ it is 
straightforward to obtain the following equation and its 
trace 
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KT-ifRRR{VRf + 2f -RJr 
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where iyRf := g''^{\J aR){'^ bR) and T := T^. Q Equa- 
tions ([3]) and Q are the basic equations we use in or- 
der to find the cosmic evolution in the model _R". We 
have employed this system of equations in the past for 
several applications, and the reader is invited to consult 
Refs. [f2, f3] for a detail discussion of this approach. Be- 
fore analyzing the cosmological situation, it is important 
to make some remarks regarding several issues that arise 
in this particular model but not in other viable f{R) 
models. In f{R) theories, one usually demands the con- 
ditions Jr > and Jrr > 0. The first one is imposed 
in order to have a positive definite effective gravitational 
constant GcS ■= Gq/ /r, while the second condition, is 
considered in order to avoid instabilities around a possi- 
ble de Sitter background [l^ . We would like to elaborate 
more about this second point. 



In we introduced the potential V{R) such that 

Vr (2/ — RfR)/3 which was relevant for tracking the 
possible de Sitter points allowed by the theory and which 
correspond to trivial solutions of Eq. Q in vacuum. 
These trivial solutions are given hy R = Ri = const. 
such that Vr{Ri) = 0, assuming fRR{Ri) ^ 0, where 
the effective cosmological constant is Aes = i?i/4. This 
explains qualitatively why f{R) theories having a de Sit- 
ter point can potentially produce an accelerated expan- 
sion when R ^ Ri and Pmatt ^ as the Universe 
evolves. Now, if fRR{Ri) = 0, one should define in- 
stead Vr = (2/ - RfR)/{3fRR), being that Jrr appears 
in the denominator of Eq. Q . Since this situation hap- 
pens generically in the model i?", we shall consider Vr 
and not Vr. Related with the stability analysis is the 
mass of scalar mode around a de Sitter point R = Ri. 

:= Vrr{Ri) = (m^ - fRRRVR/ fRR)R^, where vn? = 
Ur - RfRRUJl^lRR^Ri)], and if fRR{Ri) ^ then 
TO^ = m? since Vr{Ri) — 0. Usually when fRR{Ri) ^ 0, 
TO^ is negative if fRR^Ri) < (assuming /i?(i?i) > 
and i?i > 0), and in that case instabilities may develop 
rapidly in time [Ts'l. Thus one should consider theories 
where fRR{Ri) > 0, and in this case > if the 
critical point at i?i is a minimum of V{R) or V{R). 

Let us now focus on f{R) = XRn{R/Rn)^, where A is 
a dimensionless constant, and i?„ is a another constant 
which in general depends on the choice of n, and settles 
the built-in scale. In practice i?„ = anHg, where a„ is a 
dimensionless constant and Hq is the current Hubble pa- 
rameter. One then has Jrr = Xn{n— l)(i?/i?„)"~^i?~^. 
We shall not consider the case n — Q nor 1 because n = 1 
corresponds to general relativity (GR), for which some 
sort of dark energy or cosmological constant is required in 
order to explain the accelerated expansion, and for n = 
the theory "disappears" (i.e. it is too simple), so from 
now on we assume n 7^ 0, 1. The condition fnn > holds 
in general provided n>Iorn<— 1, assuming in both 
cases R > 0, and /rr may vanish only at i? = (we call 
this point Rq) or when i? — ^ 00 (i?oo)- We shall not con- 
sider n < because then fji = \n{R/ Rn)^^^ becomes 
negative (assuming i? > 0), and the condition Goff > 
is violated. The quantity //j also vanishes at i?o or Rooi 
depending on 71. Finally, Vr = i?,^(2 — n)/[37i(n— 1)], and 
thus V{R) = R^ {2 - n) / [9n{n -1)]+ const. E For n = 2, 



^ We assume in all the article that a subscript R stands for 8/ Sr. 



2 The model /(_R) = R-^i'^/R has a de Sitter point at _Ri = tJ?V^ 
and the mass is negative: = — /i^y^. 

^ Had we considered the potential V(R) instead oiV{K) one would 
obtain Vr(R) = = A_R„(2-n)(_R/_R„)"/3, which forn 7^ 2 and 
positive has i? = as the only stationary solution in vacuum. 
Therefore in practice V{R) and V{R) single out the same location 
for the extrema Rq and R\ which correspond to the stationary 
(trivial) vacuum solutions of Eq. ^ alluded in the main text for 
the -R" model. 
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V{R) — const., and any R ~ Ri ^ can be a de Sitter 
point, the specific value Ri depends on the initial con- 
ditions when integrating the equations. Apart from this 
"degenerate" case, Vr vanishes only at Rq. Therefore, 
for n ^ 2 the model R" does not admit de Sitter points 
and would only be able to generate an accelerated era in 
a rather transient fashion since far in the future the mat- 
ter contribution dilutes and if R reaches some equilibrium 
point it will only be at Rq which corresponds to A^g — 0. 
The mass ^ Vni^{R) = 2i?(2 - n)/[5n{n - 1)], which 
in this case is to be evaluated at R = Ri oi R = Rq (i.e. 
R = 0) vanishes identically for n = 2, regardless of the 
value of the de Sitter point Ri. Notice that rh^ — 2m^/n, 
where m was defined above. On the other hand, for n ^ 2 
the only critical point of V{R) ^ R"^ + const, is a saddle 
point at i? = {Rq), where, as mentioned before V]i{Ro) 
vanishes, and where rh^ vanishes as well regardless of the 
value of n (we assumed n ^ 0, 1). 

When a de Sitter point i?i ^ exists in vacuum Ri — 
4Aeff = 12i?4, ^ (c.f. Eq. dll) with px = Acs/n 
and in the limit p ^ 0). However, with Rq = 0, one 
is led to Acff = = H^^^^. Faraoni Q overlooked this 
fact an obtained instead — ^{Jr/ Jrr — R)r„ = (2 — 
n)i?o/[3(?i — assuming Rq ^ 0, and thus concluding 
vr? 7^ for n ^ 2. As we just argued, this conclusion is 
incorrect since the only "de Sitter" point in the i?" model 
is i?n = 0, for n^2 and therefore rr? = 0. The analysis 
in y relies on the fact that rr? ^ and requires the latter 
to be sufficiently large for the chameleon mechanism to 
ensue, in which case the author concluded n — \ with 
< (5 < 5 X 10~'^°. Again, that analysis would be valid 
if the model had a true de Sitter point at i?i 7^ for 
n ^ 2. In light of the previous discussion, we see that 
the analysis in [3| is no longer sustained nor even required 
since no matter the value of n (with n 7^ 0, 1) the scalar 
mode is massless. In reality, the chameleon requires a 
"thin shell" condition and an effective mass [l6|, both 
depending on the density of the environment, so ni? by 
its own does not suffice to analyze such mechanism. But, 
if it were the case, then the i?" model would be discarded 



Notice the missing factors of '2 and n with respect to rh? . The 
difference arises because in our definition of rh? we did not as- 
sume anything about the critical point precisely because fuR. 
might vanish there. Nonetheless, such factors are irrelevant for 
-Rq = since then m? = = fh? . 
^ In Q] the range of the scalar mode was denoted by s(n) which 
is given by s{n) = l/m ~ but since both m and m are 

zero at _Ro and at R\ for any n ^ 0,1 then s{n) — > 00, contrary 
to what was found in Q for n ^ 2, where it was assumed that 
^vac ~ ^0/12 ^ 0, denoted by in that reference. Here Hq is 
the actual cosmological constant where all forms of matter (or- 
dinary and the "geometric dark energy") are taken into account, 
while -ffvac is the Hubble expansion when the ordinary matter is 
neglected and when it is evaluated at the stationary solution of 
Eq. So in ^ no distinction is made between Hq and -ffvac- 



automatically even if 5 were within the above interval 
(with 5^0) since, the scalar mode being massless, one of 
the Post-Newtonian parameter would be 7 1/2 whose 
relative difference with 7gr = 1 is more than four orders 
of magnitude larger that the maximum value admitted 
by observations I7 - 1| < 2.3 x 10"^ '17]. 

In the next section we perform a numerical analysis of 
the full cosmological equations and show that within the 
model i?", including the case n = 2, an adequate mat- 
ter dominated era followed by a satisfactory accelerated 
expansion is very unlikely or impossible to happen. 

III. COSMOLOGY IN f{R) 

We assume a homogeneous and isotropic space-time 
described by the Friedmann-Robertson- Walker metric: 

ds^ = -dt^ + a^(t)[dr^ ^r^ {de^ ^sin^ OdLp^)] , (5) 

where we have taken the flat case k — 0. From Eqs. ^ 
and ^ we have. 



R = 


OjRR 


S/rrrR^ + 2f - JrR + kT 


,(6) 








(7) 


H = 


--ff^ - |/9 + px + 3 (prad • 


(8) 


where ' 


= d/dt and H — 


d/a, is the Hubble expansion. In 



the above equations we have included the energy density 
p associated with matter (baryons and dark matter) and 
radiation, as well as the geometric dark energy density 
Px and pressure px given explicitly by 

px ^^{UfRR-f)-3fRRHR + Kp{l-fR)], 
HJR 12 J 

(9) 

PX = -ttV (/fl^ + /) + ^RrHR -n{p- ip.^^fR) 
Sk/r [2 

(10) 

These quantities can also be obtained from a covariant 
and conserved energy-momentum tensor associated with 
the geometric modifications to GR [isj . 

Notice that the expression for the Ricci scalar com- 
puted directly from the metric ([S]) is given by i? = 



® It is important to stress that the "weak-field" , linear or Newto- 
nian limits in f{R) theories are usually studied around a maxi- 
mum or minimum of V{R). The fact that in this case the critical 
point is a saddle point indicates that a full non linear analysis 
is required around that point and that such limits are to be re- 
considered in i?" gravity. Notice that Eq. (|4]l reads explicitly 
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FIG. 1. (color online) Ricci scalar (left panel), Hubble parameter (middle panel) and the total EOS ojtot in R" gravity for several 
values of the exponent n, taking A = 1 and the constants Qn = Ru/Hq as follows: ao.9 ~ 577.85, ai.oi ~ 404.84, Q1.2 ~ 2.02, 
ai.3 ~ 1.07, Q1.5 ~ 8 X 10"*, 02 ~ 7.9 x 10"^, 0:3 ~ 2.6 x 10"®. The ACDM model is plotted for reference. The plots of the 
middle and right panels correspond to the cases of the left panel. For n < 2 the Ricci scalar and the Hubble expansion approach 
zero as z —1, while these quantities keep growing for n > 2. The model n = 2 has an effective cosmological constant and 
produces wtot = — 1 in the far future. However, it does not posses a sufficiently large matter dominated epoch with utot ~ 0, 
as one can appreciate from the right panel. None of the models decelerate and accelerate as the ACDM model. 



&{H + 2H^) which is, as one can check, compatible with 
the previous evolution equations. Therefore, one can use 
this latter instead of Eq. ([S]). The modified Friedmann 
Eq. ([7]) is used only to check the consistency and accu- 
racy of our numerical code at every time step and also 
to fix the initial data (see Ref. [l^ for the details). We 
shall not use t as independent variable but a = ln(a/ao), 
where ag is the present value of a. The corresponding 
differential equations can be found in p^ . 

The matter variables obey the conservation equation 
Pi ~ —3H {pi + Pi) for each fluid component labeled by 
i (with Pbar,DM = and p^ad = Prad/3) which inte- 
grates straightforwardly and gives rise to the usual ex- 
pression for the energy density of matter plus radiation: 
P = (Pbar+PDM)(a/ao)"^+P?ad(a/ao)"'^, where the knot- 
ted quantities indicate their values today. The X-fluid 
variables ([9]) and pO|) also satisfy a conservation equa- 
tion similar to the one above, but with an equation of 
state (EOS) ojx ■=Px/px that evolves in cosmic time. 

The different domination eras can be tracked via the 
total EOS defined by Wtot = (Prad + Px)/{p + Px) which 
using Eqs. © and ^ yields 



h UbR + /) + HbrHR - Kp 



\ UrR -!)- SIrrHR + np 



(11) 



For instance, during the radiation, matter and geometric- 
dark-energy dominated eras cjtot ^ 1/3,0,-1 respec- 
tively. Clearly such values are also correlated with the 
behavior of the dimensionless densities Q,i = npi/ (SH^) 
which satisfy the constraint f^iad + ^^matt + ^x =^ 1 where 
f^matt ■= ^^bar + ^^DM- The Capability of the i?" model 
for reproducing the correct domination eras will be as- 
sessed by the behavior of fl's and ujtot during the cosmic 
evolution relative to the ACDM model. In this regard it 
is important to remark that the X— fluid could behave as 



a matter, radiation or even as a "ghost" fluid (one with 
Px < 0) depending on the value of the exponent n, and 
therefore it could lead to an inadequate evolution history 
of the Universe. We discuss these possibilities in the next 
section. 



IV. NUMERICAL RESULTS AND DISCUSSION 

We integrate the differential equations starting at some 
redshift z = ag/a — 1, say z ^ 150, by assuming matter 
domination for all the n's in the i?" model that we ana- 
lyze. We obtain the initial conditions as described in [l^ 
and find that varying them in several ways it turns out 
impossible to recover the actual abundances of the dif- 
ferent components at present time while having an ade- 
quate accelerating phase. Here we take A = 1 but our 
conclusions do not change by choosing other (positive) 
values. This means that compared to the ACDM model, 
the Universe expands faster or slower depending on n but 
it never reproduces the correct accelerated expansion and 
matter domination eras within the same model; it repro- 
duces one or the other in the best of scenarios but not 
both. Figure [1] shows the evolution of the Hubble pa- 
rameter and the Ricci scalar from the past at z ~ 150 to 
the far future z — >■ — 1 (the current time corresponding 
to z — 0). Notice that for n = 2 the model admits a 
de Sitter solution with R Ri ^ 12^^o Universe 
evolves towards the present time. Since we have taken 
into account the matter terms, the previous equality does 
not hold exactly, but approximates very well to the ex- 
pected value, in agreement with our previous analysis of 
Secini From Figured] (right panel) we appreciate that for 
this n, the EOS oj^^^ is close to Wobs ^ —0.75, which is 
the required value to explain the current accelerated ex- 



4 



1.2 
0.9 
0.6 
0.3 


-0.3 



-1 

Log,o(1+z) 




1 

0.8 
0.6 
0.4 
0.2 




-1 1 

Log,o(1+z) 



FIG. 2. (color online) Evolution of fimatt (red solid line), Qx (blue solid line) and f^rad (green solid line) in R" gravity for 
n = 1.01 (left panel), n = 1.3 (middle panel) and n = 2 (right panel). For reference the corresponding quantities of the ACDM 
model are included in each panel (dashed lines). Notice from the left panel that Six can be negative and flmatt > 1. 



pansion of the Universe. Nevertheless, the matter epoch 
is very short as compared with the ACDM model. For 
any other value of n, a de Sitter point is never reached, 
instead i? — >■ 0, and ff— >0for0<n<2 and R and H 
grows in the future for n > 2 (c.f. equation in footnote 
6). 

The ACDM model compatible with the supernovae 
data shows that matter starts dominating for z > 0.45 
and dark energy for z < 0.45 which correspond re- 
spectively to Wtot ^ —0.5, and Wtot ^ —0.5 reaching 
wtot > -10-2 for z>5, and Wft < -0.75 for z < 0. The 
Universe starts accelerating when wtot < — 1/3 at z ~ 0.8. 
Figure [T] shows that for n ^ ric with Uc ~ 1.285 a suffi- 
ciently large matter dominated era exists with |wtot| ^ 1, 
but approaches the value OJ^^^ ~ —0.212 which is incom- 
patible with Wobs ^ —0.75. For n < Uc, there is never a 
matter domination era and wf^*^ is always far from Wobs, 
and it can even be positive. In particular, for n < 1, 
which we include for illustrative purposes as it violates 
the condition fjiji > 0, the model behaves as radiation 
dominated with utot ^1/3 and the derivatives fa and 
fnfj blow up when i? — >■ 0. Finally, for n > Uc, there 
is never a matter dominated epoch, but just a transient 
one with wtot < interpolating monotonically between 
Wtot ~ and a negative value at z = 0. Among these 
values, for n > 3 basically all the models behave iden- 
tically with wtot ~ —1.067 as z — —1. Figure [2] shows 
the corresponding evolution of the fractions fix, i^matt 
and rirad for a prototype of examples that qualitatively 
encompasses the rest of the cases, and are compared with 
the ACDM model. For n ^ ric the abundances are similar 
to ACDM, particularly at the present epoch (z = 0), but 
as we mentioned above, the model is unable to accelerate 
the Universe properly. For n > ric the matter domination 
epoch is very short (in agreement with the behavior of 
wtot in Fig. [IJ. Finally, for n < ric, i^x can even become 
negative, with fimatt possessing phases of superdomina- 
tion (i.e. fimatt > 1 in those phases) where Wtot can be- 
come positive. In particular when we take n = 1 + e with 



|e| <^ 1, the denominator in Eq. ([5]), or equivalently in 
Eq. @, becomes very small (c.f. equation in footnote 6) 
producing an important contribution on the r.h.s. of the 
differential equation for R. The cosmological evolution 
for such values of n is then rather different from general 
relativity. For instance, taking n — 1.01 we appreciate 
from Fig. [T] (right panel) that ojtot oscillates around a 
value near zero due to the oscillations of R (left panel). 
This oscillatory behavior can also be appreciated in fix 
and rimatt from Fig. [2] (left panel). The amplitude of 
these oscillations are damped so that wtot —J- in the 
present (z = 0) and future, thus the Universe does not 
accelerate. On the other hand, for e < 0, say e = —0.1 
(n = 0.9), Wtot ~ 1/3, as mentioned before, and the Uni- 
verse behaves as radiation dominated. 

Carloni et. al. [llj following ^10|], performed a cos- 
mological analysis using a dynamical system approach 
based on a first order system of equations which is dif- 
ferent from ours and which was useful for a qualitative 
description of the cosmological evolution in i?" gravity. 
In their approach they found the fixed points (stable, 
unstable or saddle) of this and other f{R) models which 
can represent the matter or the accelerated phases in the 
Universe. As argued by these authors some of these fixed 
points are different from the ones found in Q which, 
as they stressed, might change the conclusions therein. 
Nevertheless, the authors in JTi!] acknowledge that their 
analysis is only qualitative as the fixed points might not 
even be connected, and that an accurate numerical anal- 
ysis is required. This is precisely that we have performed 
here. 

In summary, the homogeneous and isotropic cosmol- 
ogy in i?" gravity seems to show a complete disagree- 
ment with what is required to explain the current fea- 
tures of the actual Universe. Since our numerical inte- 
gration was performed by including the whole mixture of 
components in the Universe, even if radiation is relatively 
small, and without any identification with a scalar-tensor 
theory in any frame whatsoever, the generic problems in 
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the i?" model seem real and are not due to any arti- 
fact concerning the ST approach or due to any inconsis- 
tency regarding a phase space analysis as objected in 
[T]| . Thus, we strongly support the conclusion that i?" is 
not a cosmologically viable candidate, unless a curvature 
k changes things dramatically and makes everything 
fit with observations. But in such occurrence, a non stan- 
dard inflationary paradigm has to be called for explaining 
the origin of cosmological perturbations. 

In [ij, we explored other f{R) models that can 
produce a successful background cosmology (i.e. without 
taking into account perturbations) but needless to say, a 



detailed scrutiny is required in all possible ambits before 
considering f{R) theories as a serious threat to general 
relativity. 
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